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Abstract 

The radiative neutrino mass model with an inert doublet scalar has been considered as 
a promising candidate which can explain neutrino masses, dark matter abundance and 
baryon number asymmetry if dark matter is identified with the lightest neutral component 
of the inert doublet. We reexamine these properties by imposing all the data of the 
neutrino oscillation, which are recently suggested by the reactor experiments. We find that 
the sufficient baryon number asymmetry seems not to be easily generated in a consistent 
way with all the data of the neutrino masses and mixing as long as the right-handed 
neutrinos are kept in TeV regions. Two possible modifications of the model are examined. 
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1 Introduction 



The standard model (SM) is now considered to be extended on the basis of several evi- 
dences clarified by recent experiments and observations, that is, the neutrino masses and 
mixing [I] , the existence of dark matter (DM) [2] and also the baryon number asymmetry 
in the Universe [3]. Although there are a lot of models which are proposed to explain 
these independently, it is not so easy to construct a model which can explain all of them 
simultaneously without causing any other phenomenological problems. If we can find such 
a model, it would give us crucial hints for the new physics beyond the SM. The study 
along this line might play a crucial role for the search of physics beyond the SM prior to 
the study based on purely theoretical motivation such as the gauge hierarchy problem. 

The radiative neutrino mass model with an inert doublet jl] could be such a promising 
candidate. It is a very simple extension of the SM by an inert doublet scalar and three 
right-handed neutrinos only. An imposed Z 2 symmetry controls the scalar potential and 
forbids the tree-level neutrino masses since its odd parity is assigned to these new fields 
and the even parity is assigned to the SM fields. It also guarantees the stability of 
the lightest field with its odd parity. Thus, the lightest neutral component of the inert 
doublet scalar [51 [6J |9] or the lightest right-handed neutrino [TOj [TT] could play a 
role of DM. This feature opens a possibility for the model such that it can explain all the 
above-mentioned three problems on the basis of closely related physics, simultaneously. 
However, the model could explain only two of the three problems if the right-handed 
neutrino is identified with DM [12]. In this case the model is required to be extended in 
some way for the explanation of all three issues [13]. On the other hand, it is noticeable 
that the above-mentioned three problems are suggested to be consistently explained by 
the new fields with TeV-scale masses as long as the lightest neutral component of the 
inert doublet scalar is identified with DM [6]. This latter case seems very interesting and 
worthy for further quantitative study since the signature of the model could be seen in 
various ongoing or future experiments. 

In this paper we reexamine this radiative neutrino mass model by fixing the parameters 
relevant to the neutrino masses and mixing on the basis of neutrino oscillation data 
including the recent results for #13 given by T2K, Double Chooz, RENO and Daya Bay [H]. 
We proceed this study by imposing the conditions required by the DM relic abundance 
and its direct detection. Based on these results, we analyze what amount of the baryon 
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number asymmetry can be generated via thermal leptogenesis. We show that the sufficient 
baryon number asymmetry seems difficult to be generated consistently for the parameters 
that are favored by the presently known phenomenological requirements. 

The remaining parts of paper are organized as follows. In section 2, we briefly review 
the scalar sector of the model and then discuss the constraints brought about by the DM 
relic abundance and the DM direct detection. Next, we fix the parameters relevant to the 
neutrino mass matrix to realize the neutrino oscillation data, which are also closely related 
to leptogenesis. In section 3, we apply them to the study of leptogenesis and estimate 
the baryon number asymmetry via the out-of-thermal-equilibrium decay of the lightest 
right-handed neutrino by solving the Boltzmann equations numerically. Conditions to 
generate the suitable baryon number asymmetry are discussed. Summary of the paper is 
given in section 4. 

2 DM abundance and neutrino masses 
2.1 The model and nature of its scalar sector 

We consider the radiative neutrino mass model with an inert doublet scalar [3] . The model 
is a very simple extension of the standard model (SM) with three right-handed neutrinos 
Ni, and a scalar doublet 77 which is called the inert doublet and assumed to have no 
vacuum expectation value. Although both iVj and r\ are supposed to have odd parity 
of an assumed Z 2 symmetry, all SM contents are assigned by its even parity. Invariant 
Yukawa couplings and scalar potential related to these new fields are summarized as 



where li is a left-handed lepton doublet and (f) is an ordinary Higgs doublet. All the 
quartic coupling constants Aj are assumed to be real, for simplicity. We also assume that 
neutrino Yukawa couplings hy are written by using the basis under which both matrices 
for Yukawa couplings of charged leptons and for masses of the right-handed neutrinos 
are real and diagonal. These neutrino Yukawa couplings are constrained by the neutrino 
oscillation data and also the lepton flavor-violating processes such as // — > ej. 



-Cy 



V 



hijNtqUj + //,y ,//.V, + -y (NiNt + N?N t ) , 



+ 




(i) 
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In the following study, we assume the mass spectrum of the right-handed neutrinos to 
satisfy 

M 1 <M 2 < M 3 , (2) 
and also the flavor structure of the neutrino Yukawa couplings to be 

foei 0? hfii foil fori Qlfoij 

h ej = hj, hpj = q 2 hj. h T j = -qahj, (3) 

where q\, 2 ,z are rea l constants. This assumption for the neutrino Yukawa couplings could 
reduce free parameters of the model substantially. Moreover, it can cause the favorable 
lepton flavor mixing as found later. We note that there remains a freedom, that is, which 
type structure represented by the suffix i and j in eq. ([3]) should be assigned to each 
right-handed neutrino. In the following part, we adopt two typical cases for it as follows, 

(i) i = 1,2, j = 3; (ii) * = 1,3, j = 2. (4) 



Now, we briefly review the scalar sector of the model [21 E]- If we take unitary gauge 
and put 4> T = (0, ((f)) + and i] T = (r] + , -^(r]R + iVi)) "where (0) = -^f-, the scalar 
potential V in eq. ([!]) can be written as 



, _ . /w,2 . /w„+„- , „2 , „2\ 
4 

1 , 2 

+ t 



X 1 h z -^X 2 (rj + rj +Vr + Vi) 



(2A 3 + 2 V / V^V + (2A+ + 2 v / V^)r/ 2 ? + (2A_ + 2 V / V^)r/ 2 J , (5) 
where we use the definition A± = A 3 + A 4 ± A 5 and 

m 2 =4A 2 <0> 2 , M 2 c = m 2 + A 3 <0> 2 , M$ R = m 2 + A + (0> 2 , M 2 , = m 2 + A„(0> 2 . (6) 
The expression of V in eq. (jHJ) shows that the assumed vacuum is stable for 



Ai, A 2 > 0, A 3 , A+, A_ > (7) 

We also require these quartic couplings to satisfy |Aj| < 4n so that the perturbativity of 
the model is guaranteed. 

Since the new doublet scalar 77 is assumed to have no vacuum expectation value, the 
Z 2 symmetry is kept as the unbroken symmetry of the model. Thus, the lightest field 
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with the odd parity of this Z2 is stable and then its thermal relic behaves as DM in the 
Universe. If it is identified with t]r here, the following condition should be satisfied 

A 4 + A 5 <0, A 5 <0; M VR <M U (8) 

These are easily found from eq. (J6j. The value of Ai might be estimated by using rrih ~ 
125 GeV, which is suggested through the recent LHC experiments. If we apply it to the 
tree-level formula in eq. (EI), we have Ai ~ 0.1. Using this value of Ai and the conditions 
given in eqs. (|7j) and (jHJ), we can roughly estimate the allowed range of A34 as 

A 3 > -1, > A 4 > -4?r, (9) 

for the sufficiently small values of [ A5 1 . The lower bound of A4 is settled by the requirement 
for the perturbativity of the model. 

The mass difference among the components of 77 is estimated as 

M m -M VR |A 5 |(0) 2 _ 5 M rh -M m |A 4 + A 5 |(0> 2 



(10) 



which could be a good approximation for the large value of m v such as 0(1) TeV§ These 
formulas show that coannihilation among the components of rj could play an important 
role in the estimation of the relic abundance of tjr [6j. 

2.2 Inert doublet dark matter 

In several articles [5], the DM abundance is found to be well explained if the lightest 
neutral component of 77 is identified with DM. In the high-mass 77 case, in particular, it is 
suggested that the relic abundance could be a suitable value if one of the quartic couplings 
|Aj| in eq. (jTJ) has magnitude of 0(1) [5J. 

The t]r relic abundance is known to be estimated as [T5] 

1 .07 x lQQGeV" 1 
J(x F )gl /2 m pl 

where the freeze-out temperature Tp(= M Vr /xf) and J(xp) are defined as 



V m h 2 ~ ' - n , (11) 



^ = ]n 0.0387n pl , eff M (a eS v) ^ J( ^ } = /~ (a^ ^ ^ 

\g*xp) 1 J Xp x 



3 Such a large value of m, is favored from the analysis of the T parameter for the precise measurements 
in the electroweak interaction [5J|S]. In that case, the model has no constraint from it. 
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The effective annihilation cross section (a e gv) and the effective degrees of freedom g e g are 
expressed by using the thermally averaged (co) annihilation cross section (a^v) and the 
r]i equilibrium number density n eq = ( ^ ) e _M, '*' T 



as 

1 ■ c q n ec| 4 ec l 

) = TLM^. ^ff = 2^^q- ( 13 ) 



9a ^ "T"i " 



The thermally averaged (co) annihilation cross section may be expanded by the ther- 
mally averaged relative velocity (v 2 ) of the annihilating fields as (<Jijv) = ciij + bij(v 2 ). 
Since (v 2 ) <§C 1 is satisfied for cold DM candidates and then gives the dominant role for 
determining the relic abundance of 7/^, we take account of it only neglecting the 6^- con- 
tribution in this analysis. In the present model, the corresponding cross section is caused 
by the weak gauge interactions and also the quartic couplings \. It is approximately 
calculated as [6] 

+ sM^ {Nl3 + Nu + N23 + Nu) 

+ [(A+ + A 2 + 2X 2 3 )(N 11 + N 22 ) + (A+ - A_) 2 (7V 33 + iV 44 + N 12 ) 

+ {(A + - A 3 ) 2 + (A_ - A3) 2 } (JV a3 + N u + N 23 + N 24 ) 

+ {(A + + A_) 2 + 4A 2 }iV34], (14) 



where Ny is defined as 



1 n eq n eq 1 f M M \ 3/2 

— „2 eq eq — o I Ji >r9 I c ^ 



M 2 
m 



M Vi + M Vj - 2M V1 
T 



(15) 



Using these formulas, we examine the condition on the relevant parameters of the model 
to realize the relic abundance ^l VR h 2 = 0.11 which is required from the WMAP data [2]. 
In Fig. 1, we plot Q VR h 2 as a function of A4 for A5 = — 10~ 5 and some typical values of 
A3. We should note that the allowed region of A4 is restricted by eq. ([9]). This figure 
shows that the required value for fi VR h 2 could be obtained for a wide range value of m v if 
IA3 + A4I has a value of 0(1). Since we consider the high- mass region such as m v ^> ((f)), the 
coannihilation among the components of rj could be effective to reduce the relic abundance 

of T) R - 



We may use the notation such as rj 2 , 773, 774) = (^i?, 777, f? ) for convenience in the following 
discussion. 
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Fig. 1 Sl nR h 2 as a function of A4 for the negative value of A3 such as —1, — 0.5, — 0.1 in the left panel 
and for the positive value of A3 such as 1, 0.5, 0.1 in the right panel. The value of m, is fixed to 1000 
and 1500 GeV in both cases. 

The above analysis shows that the DM relic abundance gives only a weak condition 
on some of the quartic couplings but no conditions on the neutrino Yukawa couplings. 
This is completely different from the case in which the lightest iVj is identified with DM 
P21 IHJ EE] • On the other hand, the direct search of DM could give a severe constraint 
on the value of A5, which plays a crucial role in this radiative neutrino mass generation. 
Elastic scatterings between r] R and nucleus could be mediated by the Higgs exchange at 
tree level and also by the gauge boson exchange at one-loop level. However, their effects 
are much smaller than the present upper bounds of the sensitivity for the direct detection. 
Thus, we can neglect their effects in this discussion. On the other hand, inelastic scattering 
of rjji with nucleus mediated by the Z° exchange could bring about an important effect to 
the direct search experiments [8j [9] , since the masses of rjR and rji are almost degenerate 
for small values of | A5 1 as found from eq. (TIP]) . 

If we note that the interaction of r] R relevant to this process is given by 



£ = griRd^rjiZ 11 - gr)id fJ ,r] R Z IJ ', 



(16) 



it is found that the inelastic nucleus-DM scattering can occur for the DM with velocity 
larger than a minimum value given by [T6] 



1 



(rn I1 E IL + 5 



p — r- \ / ( 17 ) 

where 8 is the mass difference between rjn and rjj defined in eq. ljTOl) . Er is the nucleus 
recoil energy, and m N and /i^r are the mass of the target nucleus and the reduced mass 



of the nucleus-DM system. Thus, the mass difference 5 is constrained by the fact that 

n 

no DM signal is found in the direct DM search yet [T71 dB]Jj This condition might be 
estimated as 5 ^150 keV [9]. Since S is related to A5 through eq. fllO]) , this constrains the 
allowed value of |A 5 | such aq^ 

We take account of this constraint in the following analysis of the neutrino masses and 
the baryon number asymmetry. 



2.3 Neutrino masses and mixing 

Neutrino masses are generated through one-loop diagrams with the contribution of new 
Z 2 odd fields. They can be expressed as [U HO] 



k=l L '/ ft \ '/ k ' k=l 

where = + (A 3 + A 4 ) (0) 2 . Since we consider the high-mass region such as ^> ((f)), 
the mass difference among rji caused by nonzero A4 is negligible in the neutrino mass 
analysis. Thus, we treat their masses as M v . Both neutrino masses and mixing are 
determined by the couplings A5 and hik, the right-handed neutrino masses Mj's and the 
inert doublet mass M n . Here we note that the neutrino Yukawa couplings could take 
rather large values even for the light right-handed neutrinos with masses of 0(1) TeV as 
long as IA5I takes a small value in the range given by eq. (fl8l) . This freedom is crucial 
when we consider leptogenesis in this model as seen in the next section. 

Now we have a lot of information on the feature of lepton flavor mixing on the basis 
of the neutrino oscillation data including the recent results for # 13 [19]. We can use it 
to restrict the neutrino Yukawa couplings. The flavor structure of the neutrino Yukawa 



5 The DAMA data have been suggested to be explained by the DM inelastic scattering [8, 9 . However, 

we do not consider it here. 

6 We should note that the bound of 6 largely depends on the DM velocity in the neighborhood of the 

Earth. If we take 5 ^1 MeV, this inelastic scattering effect can be completely neglected even for the 

maximally estimated DM velocity. In that case, the lower bound of |As| becomes one order of magnitude 

larger. 
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1 qi 
V Qi Qi ) 



{h\A x + hjK) + 



(20) 



couplings assumed in eq. ([3]) makes the neutrino mass matrix take a simple form such as 
/ \ / 1 q 2 -q 3 \ 

q 2 q\ -q2qz 

\ -<?3 -<72<?3 

where i, j should be understood to stand for (i) i = 2, j = 3 and (ii) i = 3, j = 2 following 
eq. (jl]). If we put (71,2,3 = 1 in both cases, the PMNS mixing matrix is easily found to 
have a tribimaximal form 

\ / 1 \ 



U 



PMNS 



/ 2 1 

-1 1 1 
VG V3 a/2 







(21) 



/ 



e 

where Majorana phases 0:1,2 are determined by the phases hi and A5. If we put ^ = arg(/i 
and ip\ 5 = arg(As), they are expressed as 



«2 



1 

- tan 
2 



I'll 


2 A 1 sin(2(^ 1 + ip Xs ) + 


/l2 


2 A 2 sin(2v9 2 + ip X5 ) 


/tl 2 Ai COs(2(y9i + V?A 5 ) + |^-2| 2 A 2 COs(2v?2 + V^Ag) 



■ (22) 



In this case, one of mass eigenvalues is zero. Thus, if \hi\ is assumed to take a sufficiently 
small value compared with others, we find that the mass eigenvalues should satisfy 

\/Am 2 s 



2 A _yA§£ : 



\hA A4 



l sol 



(23) 



where Am 2 tm and Am 2 ol stand for the squared mass differences required by the neutrino 
oscillation analysis for both atmospheric and solar neutrinos [TJ [19] 

Using the formulas (T2~T]) and ff23|) . we can examine whether parameters obtained in 
the previous part could be consistent with the neutrino oscillation data. This gives a 
useful starting point for the analysis. However, if we take account of the fact that #13 is 
found to have a nonzero value now, we cannot use them in the analysis directly. Here, we 
numerically diagonalize the mass matrix ( 120]) and impose both the neutrino oscillation 
data and the constraint on |As| given in eq. (llSp to restrict the neutrino Yukawa couplings. 
By fixing qx to typical values, in Fig. 2, we plot contours in the (52,53) plane which 
correspond to 2a bounds of the neutrino oscillation parameters given in [20J . This figure 
shows the mass matrix (|20|) can explain all the neutrino oscillation data consistently at 



7 We confine our study to the normal hierarchy here. 
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-0 .5 . 5 1 1. 5 2 -0 .5 . 5 1 1. 5 2 

q2 q2 

Fig. 2 Regions in the (<72, 93) plane allowed by the neutrino oscillation data for qi = 0.85 (left panel) 
and qi — 1 (right panel), which are contained in the circle drawn by the dotted line. Relevant parameters 
are fixed to the ones shown in Table 1. Each contour in both panels represents 2<r boundary values 
of neutrino oscillation parameters Am§ 2 (thick red solid and dashed lines), |Am 2 2 | (thin red solid and 
dashed lines), sin 2 2# 2 3 (greed solid and dashed lines), sin 2 2#i2 (blue solid and dashed lines) which are 
given in Ref.(3D]. The 90% CL value of sin 2 2#i3 given in Ref.JT^ is also plotted as a reference (black 
solid and dashed lines). 

the regions in the (52,93) plane, which are the region including (—0.27,2.4) in the left 
panel (qi = 0.85) and the region including (0.05,2.1) in the right panel (qi = 1). These 
regions in two panels are obtained for each set of parameters listed in Table 1. We also 
give the predicted values for sin 2 2#i3 in each case of the same Table. These examples 
show that neutrino Yukawa couplings of O(10~ 3 ) can explain the neutrino oscillation data 
for IA5I = O(10~ 5 ) and the right-handed neutrinos with the mass of 0(1) TeV. 

Lepton flavor-violating processes such as \i — > e^y are also induced through one-loop 
diagrams which have rj and iV, in the internal lines pU]. Their present experimental 
bounds could impose severe constraints on the model depending on the values of neutrino 
Yukawa couplings h^. However, since the small \hi\ of O(10 -3 ) can realize the appropriate 
values for neutrino masses even for the TeV-scale values of Mj and M v as discussed above, 
the new contributions to the lepton flavor-violating processes are sufficiently suppressed 
such as Br(/i — > ej) = O(10~ 18 ) and Br(r — > fij) = O(10~ 14 ). These values show that 
the lepton flavor- violating processes bring about no substantial constraints on the model. 
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Qi 


M v 


Mi 


M 2 


M 3 


10 3 \h 2 \ 


I0 3 \h 3 \ 


sin 2 26 13 


max £2,3 




Y B 




(ia) 


0.85 


1 


2 


6 


10 


3.41 


1.50 


0.085 


1.1 • 10~ 7 


2.7 


■ 10" 


-12 


(ib) 


0.85 


1 


2 


20 


200 


4.62 


4.16 


0.085 


6.1 • 10~ 8 


1.6 


• 10" 


-12 


(ic) 


1 


1 


2 


6 


10 


3.41 


1.50 


0.053 


9.8 • 10~ 8 


2.8 


• 10- 


-12 


(iia) 


0.85 


1 


2 


6 


10 


1.34 


3.81 


0.085 


1.7 - 10- 8 


2.7 


• 10- 


-13 


(iib) 


0.85 


1 


2 


20 


200 


1.82 


10.6 


0.085 


9.5 • 10~ 9 


8.3 


• 10" 


-13 



Table 1 The predicted value of sin 2 29\z and Yb for the model parameters that can satisfy the neutrino 
oscillation data. Cases (i) and (ii) correspond to the ones defined in eq. ©. In all cases, |As| and \h\\ are 
fixed to 10 -5 and 3 • 10 -8 , respectively. The value of sin 2 2^13 is evaluated at (92,93) = (—0.27,2.4) for 
qi = 0.85 and (0.05, 2.1) for q\ — 1, where all other neutrino oscillation data are satisfied. A TeV unit is 
used as the mass scale. 

We should note that the freedom of A5 in this mass generation scheme makes it possible. 



3 Baryon number asymmetry 

3.1 Leptogenesis via the decay of the TeV-scale right-handed 
neutrino 

We consider the thermal leptogenesis |21j in this model with the mass spectrum given 
in eqs. (|2J) and ©. In this case, the lepton number asymmetry is expected to be gener- 
ated through the out-of-thermal-equilibrium decay of the right-handed neutrino N\. The 
dominant contribution to the CP asymmetry e in this decay is brought about by the in- 
terference between the tree diagram and the one-loop vertex diagram as usual. However, 
we should note that the 77 mass is not negligible compared with the one of N\ in this 
model. Taking account of this feature, e can be calculated as [22] 



167T 

= £ 2 sin2(^i 



3 + 1 f 1 _ M v 

4 ~ 4 I Iff 



E 

i=2,3 



Im 



E 



k=e,/i,T 



hkih* ki 



12k=e,fi,T hklh* kl 



G 



Ml" Ml 



ip 2 ) + e 3 sin2(^i - (p 3 ), 



(24) 
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where G(x, y) is defined by 

G(x, y) = ^F{x, 0) + ±F(x, y) + \ (1 - yf [F(x, 0) + F(x, y)} , (25) 
and F(x, y) is represented as 

1 — y + x 



F(x,y) = yfx 



1 — y — (1 + x) In 



a; 



(26) 



If we use the flavor structure of neutrino Yukawa couplings (i) and (ii) given in eq. 
e 2 ,3 are expressed for each 



m e -C(l + a>)\h?Gl M * M A e - ° {q2 - ( M l M v 

(1) 1 - C ^ + ^\ G l Mf' MfJ' £3 ~ TTif WW 

C{q2-qiqz) 2 \h 2 \ 2 n (M 2 2 M*\ (Mi M 2 



£2 ~ iTgf G [w^)' 3 ~ { Ql)l sl \Wm!) 



,(27) 



where C = I6ir 



4 ~ 4 \^ Mf J 



The decay of Ni should be out of equilibrium so that the lepton number asymmetry is 
generated through it. If we express the Hubble parameter and the decay width of N\ by 
H and rjy respectively, this condition is given as H > T£ at T ~ Mi where the lepton 
number asymmetry is considered to be dominantly generated. Since r£ is expressed 
as = + Qi) 2 Mi ^1 — , we find that the Yukawa coupling \hi\ should be 

sufficiently small such as 

|/, 1 |<2xlO- 8 (l + g 1 2 )- 1/2 ( T ^ V ) 1/2 - (28) 

We note that this constraint could be weaker since both M v and the Boltzmann suppres- 
sion factor are neglected in this estimation. As found in the numerical calculation, \hi\ 
can be somewhat larger than this bound. 

The generated lepton number asymmetry could be washed out by both the lepton 
number violating 2-2 scattering such as 7777 — > £ a £p and r]£ a — > fyl^ and also the inverse 
decay of N\. If the relevant Yukawa couplings are much smaller than 0(1), these processes 
are expected to decouple before the temperature T of the thermal plasma decreases to T ~ 
Mi. In order to study this quantitatively, we numerically solve the coupled Boltzmann 
equations for the number density of Ni and the lepton number asymmetry which are 
expressed by and hl here, respectively. The Boltzmann equations for these quantities 
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^+^ 3 ')\, (29) 



are written as [23] 

dY Nl 
dz 

dY L _ " j / ) \, , \ . A, -U / 12. . ( !.:-!) 

dz 

where z = 4^f and H(Mi) = 1.66(7*^ 2 ^^-. and Yl are defined as = and 
Y L = — by using the entropy density s. Their equilibrium values are expressed as 
Y^{z) = 2n^g r z2 -^ 2 ^ an< ^ ^ q = wnere i s the number of relativistic degrees of 
freedom and K 2 (z) is the modified Bessel function of the second kind. In these equations 
we omit terms whose contributions are considered to be negligible compared with others. 
The formulas of the relevant reaction density 7 contained in these equations are given in 
the Appendix. If we use the relation B = ^(B — L) which is derived with the chemical 
equilibrium condition in this model, the baryon number asymmetry Yb(= — ) in the 
present Universe is found to be estimated as 

Yb = -^Y L (z EW ) (30) 

by using the solution Y L {z) of the coupled equations in eq. Here zew is related to 

the sphaleron decoupling temperature Tew as ^ew = 

The solutions of eq. (|29|) are shown in the upper panels of the left and the middle 
columns in Fig. 3 for cases (ia) and (ib) in Table 1. The generated lepton number asym- 
metry \Yl\ is smaller than the value required for the explanation of the baryon number 
asymmetry at least by one order of magnitude,^) In case (ia), the CP asymmetry param- 
eter |e| takes rather large value such as 1.1 ■ 10~ 7 . This suggests that the washout of the 
generated lepton number asymmetry is effective. In the figure of the right column we use 
the same parameters as the one of case (ia) except for \h\\ which is fixed to the larger 
value 5 • 10 -7 . From this figure, we can see the role of this coupling which is discussed 
above. Although this change does not affect the values of [ /t2,3 1 which explain the neu- 
trino oscillation data, we expect that the deviation of the number density of N\ from the 
equilibrium value becomes smaller than other cases with the smaller value of \h\\. This 
is shown in the figures. 



8 The parameters used in case (ia) are almost equivalent to the one which is presented as the promising 
one for the generation of the sufficient baryon number asymmetry in [B] . 
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Fig. 3 The upper panels show the evolution of Yl, Ym 1 and Ajvi = \Y^ 1 — Y^\. The lower panels show 
the reaction rates of the processes that have crucial effects for the leptogenesis in this model. We use 
the parameters shown as (ia) and (ib) in Table 1 for the left and middle panels. In the right panel the 
same parameters as (ia) are used except for \h\\, which is fixed to \h\\ = 5 ■ 1CP 7 in this case. The black 
dotted line represents the required value for \Y^\. 

In the lower panel we plot the behavior of the relevant reaction rates for each case. 
These processes are crucial for the leptogenesis in this model. The figures show that the 
lepton number- violating scatterings induced by the s-channel iVj exchange are kept in the 
thermal equilibrium until rather late period and the large part of the generated lepton 
number asymmetry is washed out. This situation is common for all cases in Table 1. 
We give the predicted value of Yb for each case in the last column of Table 1. These 
examples show that the sufficient amount of baryon number asymmetry seems difficult to 
be generated through the thermal leptogenesis in the present neutrino mass generation 
scheme at least as long as we impose the full neutrino oscillation data and the DM direct 
search constraint. 

In order to confirm this statement in case (ia), we plot the generated baryon number 
asymmetry Yb for various values of |As| in Fig. 4. The maximum value of Yb is found to 
be realized at a certain value of | A5 1 and it moves to the smaller region for the smaller 
M\. This may be explained as follows. At the region with a larger \\^\ value, the neutrino 
Yukawa couplings have smaller values to give a smaller value for |e|. On the other hand, 
at the region with a smaller \X 5 \ value, the neutrino Yukawa couplings have larger values 
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Fig. 4 The dependence of the generated baryon number asymmetry Yb on | As | in case (ia) with the 
different value of M\. The value of M\ is fixed to 1.3 TeV, 2 TeV and 3 TeV in each case. 



to bring about the large washout. For the Mi = 1.3 TeV case, the requited Yb can be 
obtained at a rather small value such as < 10 -6 due to this nature. However, it is 
excluded by the direct search experiments as shown in eq. (|T8|) FI Although we do not 
search the whole parameter space, we could say that the above-mentioned result does not 
change so easily for the values of |Ab| which satisfy the condition (TTgj) . This is suggested 
by the fact that the value of Yb becomes smaller for the larger |As| which makes the 
neutrino Yukawa couplings smaller under the constraints of the neutrino oscillation data. 

Finally, we give some comments on the leptogenesis in the case where the right-handed 
neutrinos have large masses comparable to the ones in the ordinary seesaw caseo In this 
case, we find that the DM relic abundance and the neutrino masses and mixing could be 
explained consistently by setting \X 5 \ and the neutrino Yukawa couplings appropriately. 
The right-handed neutrino masses could be smaller than the ones in the ordinary tree- 
level seesaw case for the same neutrino Yukawa couplings since the neutrino masses are 
generated through the one- loop effect. In the ordinary type I seesaw scenario with the 
hierarchical right-handed neutrino masses, there is an upper bound for the CP asymmetry 
which is known as the Davidson- Ibarra (DI) bound [25] and may be written as |edi| = 
sir ~~7^y^ ^- I n ^ ne P resen t case with the assumed flavour structure, the CP asymmetry 
|e| is related to the DI bound as \e\ = l g Di| 2,\n{M 2 /M ) "fxff' (iff) ' This snows that the DI 

9 It is useful to note that the smaller M v brings about the smaller value for the lower bound of |As|. 

However, its effect is only a change of the factor in case of the high-mass rj. 

10 This is considered in [7]- However, the neutrino oscillation data are not imposed in a quantitative 

way there. 
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bound could be evaded depending on the value of [ A5 [ . However, this feature does not 
mean that the model causes more efficient leptogenesis than the ordinary seesaw model. 
Since the smaller | A5 1 requires the larger neutrino Yukawa couplings under the constraints 
of the neutrino oscillation data, the washout of the lepton number asymmetry is expected 
to become large. 

We examine this aspect under the assumed lepton flavor structure by imposing all 
the neutrino oscillation data quantitatively. We assume the hierarchical right-handed 
neutrino mass spectrum and fix the parameters as follows, 

\h 1 \ = 10'\ M! = 10 Q GeV, M 2 = 10 Q+1 GeV, M 3 = 10 Q+2 GeV, (31) 

where \hi\ is determined by taking account of the condition (|28|) . We find that the CP 
asymmetry \e\ can be written for these parameters as |e| ~ [£di| (a_2)|A s | where the DI 
bound I^diI is given as |£di| < 1-9 • 10 _16+a . This relation shows that the CP asymmetry 
could escape the DI bound without causing the contradiction with the neutrino oscillation 
data if \\$\ and Mi take suitable values. In order to fix the values of neutrino Yukawa 
couplings, we impose the neutrino oscillation data in the same way as in the previous 
examples with qi = 0.85, q 2 = —0.27, and q% = 2.4. For such neutrino Yukawa couplings, 
we obtain sin 2 2# 13 = 0.085 independently on the value of \Xg\. The baryon number 
asymmetry Yg obtained at z = 20 through the analysis of the Boltzmann equations is 
plotted for some typical values of |As| and Mi in the left panel of Fig. 5. Since this z 
is much smaller than z^w — Ys(oo) could be much smaller than the plotted value 
if the washout effects do not decouple. However, we can confirm that the plotted Yb is 
recognized as Yb (00), at least for |As| > |A!f ax |, where |A™ ax | gives the maximum value of 
Yb ■ In the region of | A5 1 < |A™ ax |, the neutrino Yukawa couplings become large enough to 
continue reducing the generated lepton number asymmetry through the lepton number- 
violating scatterings which do not decouple at z = 20 completely. Since \h 2 ^\ is required 
to be larger for the smaller values of | A5 1 from the neutrino oscillation data, the large part 
of the generated lepton number asymmetry is considered to be washed out effectively 
although a larger |e| value is expected. The figure shows that the required value of Y B 
is generated for Mi > 10 8 GeV, which is somewhat smaller than the one required in the 
ordinary seesaw casein 



n The required baryon number asymmetry could be generated even in the case with Mi = O(10 8 ) GeV 
if special texture for the neutrino mass matrix is assumed even in the ordinary seesaw case [55] . 
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Fig. 5 The baryon number asymmetry Yb in cases of the heavy right-handed neutrinos (the left panel) 
and the degenerate right-handed neutrinos (the right panel). In the left panel, a GeV unit is used for the 
mass scale. In the right panel, M\ = 2 TeV is assumed. 

3.2 Improvement by suppressing the washout 

In the previous part we found that it is difficult to generate the sufficient baryon number 
asymmetry in consistent with all the neutrino oscillation data and the DM direct search 
as long as the mass of the lightest right-handed neutrino is assumed in a TeV range. 
This result is considered to be brought about by the large washout effect of the generated 
lepton number asymmetry. Here, we consider a possible improvement of this situation 
by making the neutrino Yukawa couplings small enough to suppress the washout. In this 
improvement, the CP asymmetry |e| should be kept to a suitable value such as O(10 -7 ) 
or more, simultaneously. Resonant leptogenesis can realize it. 

As such an example, we suppose that the right-handed neutrino masses Mi and M<i 
are nearly degenerate in case (ia). M\ is fixed to 2 TeV and M 2 is replaced with M 2 = 
(1 + A) Mi. In this case, we can make the neutrino Yukawa couplings much smaller 
than the ones in case (ia) by assuming a larger value for |As|. For instance, if we put 
| A5 1 = 10~ 3 , all the neutrino oscillation data can be satisfied for \h,2\ = 3.1 • 10 -4 and 
\h 3 \ = 1.5 • 10~ 4 . Although the smaller neutrino Yukawa couplings tend to make the 
CP asymmetry \e\ smaller, we can enhance |e| by supposing the degenerate right-handed 
neutrino masses. On the other hand, in that case the washout effect could be suppressed 
sufficiently. In the usual resonant leptogenesis with the right-handed neutrino masses at 
TeV regions [22 | 127 ] 128] . the generation of the sufficient lepton number asymmetry requires 
rather strict mass degeneracy. We examine whether this situation can be changed in this 
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neutrino mass generation scenario. 

In the nearly degenerate right-handed neutrino case, the dominant contribution to the 
CP asymmetry in the decay of Ni comes from the interference between the tree and the 
self-energy diagramso In the present case, the CP asymmetry \e\ can be expressed as 

(Ml-Ml)M l Y 2 . 2Af 2 . 

sin2(^ - <p 2 ) ~ - — — — =-sin2(^i - (pi), (32) 



{Ml - Mf) 2 + Mfr| vri 4A 2 + f 2 

1M!/1 _l„2N f, MV ' 1 



where T 2 = + (l — • F° r case (ia) with sin2(<^i — tp 2 ) = 0(1), |e| can have 

a value of 0(1) for A ~ 5 ■ 10~ 7 . Thus, we can expect that the mass degeneracy required 
to bring about the sufficient baryon number asymmetry could be much milder than the 
usually assumed value A ^ 1CT 8 in the TeV-scale resonant leptogenesis [28J. 

We use the formula for e given above and estimate Yg by solving the Boltzmann 
equations in eq. (129]) for some typical values of A by varying the value of | A5 1 . In the right 
panel of Fig. 5, we plot the numerical results of the obtained baryon number asymmetry. 
The required baryon number asymmetry is found to be generated for the right-handed 
neutrinos with the mass degeneracy A = O(10~ 4 ) for each value of | A5 1 . This degeneracy 
is much milder in comparison with the ordinary resonant leptogenesis at TeV scales. The 
result is caused by the nature of the model such that the neutrino Yukawa couplings 
can take sufficiently small values to suppress the washout effect keeping the neutrino 
masses in the appropriate range for the explanation of the neutrino oscillation data. The 
freedom of A5 makes it possible. The present neutrino mass generation scheme can give 
the consistent explanation for the three phenomenological problems in the SM if only 
the rather mild mass degeneracy between two light right-handed neutrinos is assumed 
without large extension of the model. 



12 Since the Yukawa couplings of 7V2 ; 3 are required to be much larger than = 3 • 10 -8 for N±, A^2,3 
are considered to be in the thermal equilibrium. Thus, we can not expect a substantial contribution to 
the lcpton number asymmetry from the decays of ^2,3. 
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4 Summary 



The inert doublet model extended with the right-handed neutrinos is a simple and inter- 
esting framework for both neutrino masses and dark matter. In this scenario, the lightest 
right-handed neutrino or the lightest neutral component of the inert doublet can be a 
dark matter candidate. Since the neutrino Yukawa couplings should be 0(1) to reduce 
the relic abundance of dark matter in the former case, the thermal leptogenesis is difficult 
due to the strong washout effect. On the other hand, the neutrino Yukawa couplings can 
be irrelevant to the dark matter abundance in the latter case. Thus, the latter scenario 
has been considered to give the consistent explanation for the origin of the small neutrino 
masses, the existence of dark matter and the baryon number asymmetry in the Universe. 

In this paper, we reexamined the possibility for the simultaneous explanation of these 
three problems in the latter case. We took account of the quantitative explanation for 
all the neutrino oscillation data including the recent results for sin 2 2# 13 in the analysis. 
The results of our study suggest that the sufficient amount of baryon number asymmetry 
seems not to be generated in a consistent way with the full neutrino oscillation data. 
The neutrino Yukawa couplings could be large enough to enhance the CP asymmetry 
even for the 0(1) TeV right-handed neutrinos in the consistent way with the neutrino 
oscillation data by using the freedom of A5. However, the same Yukawa couplings induce 
the lepton number-violating scattering processes to wash out the generated lepton number 
asymmetry. We checked this point through the numerical study. Although we do not 
study the whole parameter space, this feature seems to be rather general and then it seems 
not so easy to find the parameters to escape this situation without any modification. 

We also examined the same problem in the modified situation such as the case with 
the heavy right-handed neutrinos like the ordinary type I seesaw and also the case with 
the degenerate light right-handed neutrinos. In the former case, we found that the right- 
handed neutrino could be somewhat lighter than the ordinary seesaw one but needs to be 
heavy enough to the similar level to it. In the latter case, even if the right-handed neutrino 
masses are in a TeV range, the resonant effect can enhance the CP asymmetry to generate 
the sufficient amount of baryon number asymmetry even for the small neutrino Yukawa 
couplings which can suppress the washout sufficiently. An interesting point is that the 
required mass degeneracy in this case is much milder than the one in the ordinary seesaw 
case. This possibility is brought about by the characteristic feature in this inert doublet 
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model with the radiative neutrino mass generation. More complete study of the parameter 
regions which are not searched in this paper will be presented in future publication. 
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Appendix 



In this Appendix, we give the formulas of the reaction density contributing to the Boltz- 
mann equations for the number density of Aq and the lepton number asymmetry. For 
the processes relevant to their evolution, we could refer to the reaction density given in 
|24j . In the present model, however, interaction terms of 77 and Aq are restricted by the 
Z 2 symmetry and also the masses of 77 and Ni take the similar order values, which cause 
large difference from the ordinary seesaw leptogenesis. Thus, we need to modify these 
formulas by taking account of the features of the present modelj^ 

In order to give the expression for the reaction density of the relevant processes, we 
introduce dimensionless variables 



x 



M 2) a ° MY ~ v Mf 
where s is the squared center of mass energy. The reaction density for the decay of Aq 
can be expressed as 



M 2 
77 



(33) 



7? 



'1 + Qi)\hi\\^ „ Mz) 



-M*(l - a v ) 



An 3 ' ' z 

where K\(z) is the modified Bessel function of the first kind. 

The reaction density for the scattering processes is expressed as 

T 



(34) 



T f°° ( v/s\ 

l(ab -> ij) = -—^ / ds &(s)y/sK! ( — J 



(35) 



where s min = max[(m„ + m b ) 2 , (m, : + nij) 2 } and a(s) is the reduced cross section. In order 
to give the expression for the reaction density of the processes relevant to eq. ( 129|) . we 
define the following quantities for convenience: 



1 



DAx) 



\x - {-s/Oi + 

X 

In 



+ afci ' 



dj + 2a v + 



L 'ij 



In 



x — a,i — cij + 2d v 

y/x(x — d 




dj - 2d v ) + \J\ij(x 



Ad r 



t/x(x — di — dj — 2d n ) — a/ \ij(x — Ad r] ) 



(36) 



13 Although the modified ones are given in Appendix of |12j . the mass spectrum assumed there is 
different from the present one. The following formulas are arranged to applicable to the scenario in this 
paper. 
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As the lepton number violating scattering processes induced through the N{ exchange, 
we have 



1 (x- a v ) 2 



2tt 



3 

E 



x 



00 I 00 Qji 



X 



[x — a r , 



DAx) 



In 



2x (x — a v ) 2 
+ Di(x) + 2D i (x) 2 

x(x + a,-, — 2a r 



i>j 



+ 
+ 



X 



x 2 

(x - a n ) 2 



x 

2(x + cii 



x 



x 
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00 Qj'i OC CL f 



+ 



(x — CLi)(x — Cij) 



2a r 



CL j Q j % 



x + a, — 2a r 



In 



x(x + cii — 2a v ) 



• I ■ (1/ 



2(x + cij — 2ci v ) x + cij — 2a v 



Q"i Q j 
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In — - — - 



OCQ/rj CL, 



n 



(37) 



for l a ry — > £pr) and also 



a 



(13), 
N 1 



X 



1 

2^ 
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i=i 

CI; 



x + 2<2j — 2a r 



In 



a 2 (x 2 — \.xa^) x l 2 

Q^ 00 | (^Q^ Qy^ 

x + (x 2 - 4xa,,) 1/2 + 2a» 



2a T 



2 - 4xa n y/ 2 + 2cii - 2a r 



+ J2 ^[( hh %]y^j ln 



x + a, + a,- 



2a r , 



x + (x 2 — Axa^) 1 ! 2 + a; + 



2a r 



4x6^) x / 2 + a, + a, — 2( 



(3? 



for l a £p — )■ 7777. Here we note that cross terms has no contribution if the maximum CP 
phases are assumed in the way as sin2(<yj 2 ,3 — v?i) = 1 with <pi = 0. We adopt this 
possibility in the numerical analysis, for simplicity. Since another type of lepton number 
violating processes A^iVj —> t a tp induced by the r\ exchange have additional suppression 
due to a small |As|, we can neglect them safelyj^ 

As the lepton number conserving scattering processes which contribute to determine 
the number density of Ni, we have 

\ ,/X7/ 

1 + 



AT. AT. (x) 
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47T 



+ 



J2 \(hh%(hrf) 
aj + a,- — 2a 
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x 



a v )(a,j - a v ) 



x 



-Lij - Re[(hh% 



(cii — CL v ){aj — a v ) + xa v 
2 ^yci{CijL ij 



CLj I 2ci^ 



(39) 



14 We should note that |Ag| might not have a small value in the case with heavy right-handed neutrinos, 
which is discussed in section 3. In this case, these processes could give large contribution to the washout 
of the generated lepton number asymmetry. 



22 



for NiNj ->■ 



a 



(3) 



l a £p which are induced through the r/ exchange and also 

, , 1 
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4a v (ai — ctj) 



+ l 



Re[(^ t ) 
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(a; 2 — 4xa v y/ 2 (x — a. t — aj — 2a r . 



(40) 



for N{Nj — >■ 777^ which are induced through the £ Q exchange. The cross terms in these 
reduced cross sections are neglected because of the same reasoning as eqs (13"?]) and (}38l) . 

In order to see the behavior of these relevant processes such as the decoupling time, 
we may estimate the ratio of the reaction rate to the Hubble rate ^ as a function of z (see 
Fig. 3). The thermally averaged reaction rate V is related to the above discussed reaction 
densities through 



for the decay of N\ and also 



,(2,13) 
N 



piVi 
1 D 



(2,13) 

Tn 



cq 



(41) 
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eq ' 



F (2,3) 
1 NiN! 
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eq 



(42) 



for the 2-2 scattering processes given in eqs. fl37j) . fl38|) and eqs. (1391 . (T40]), respectively. 
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